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ABSTRACT 

We propose and test a scenario for the assembly and evolution of luminous matter in 
galaxies which substantially differs from that adopted by other semianalytic models. As 
for the dark matter (DM) , we follow the detailed evolution of halos within the canonical 
ACDM cosmology using standard Montecarlo methods. However, when overlaying 
prescriptions for baryon evolution, we take into account an effect pointed out in the 
past few years by a number of studies mostly based on intensive iV-body simulations, 
namely that typical halo growth occurs in two phases: an early, fast collapse phase 
featuring several major merger events, followed by a late, quiescent accretion onto 
the halo outskirts. We propose that the two modes of halo growth drive two distinct 
modes for the evolution of baryonic matter, favoring the development of the spheroidal 
and disc components of galaxies, respectively. We test this idea using the semianalytic 
technique. Our galaxy formation model envisages an early coevolution of spheroids and 
the central supermassive black holes, already tested in our previous works, followed 
by a relatively quiescent growth of discs around the preformed spheroids. In this 
exploratory study, we couple our model to the spectrophotometric code GRASIL, 
and compare our results on several properties of the local galaxy population with 
observations, finding an encouraging agreement. 
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1 INTRODUCTION 

A fundamental issue when modeling the evolution of galax- 
ies in a cosmological context is that the majority of the pro- 
cesses driving baryonic evolution (such as star formation, 
various feedback mechanisms, accretion onto supermassive 
black holes [BHs]) operate or originate on scales well below 
the resolution of any feasible simulation in a cosmic box. 
Moreover, these processes are highly nonlinear, poorly un- 
derstood from a physical point of view, and approximated by 
means of simplified, often phenomenological, and thus un- 
certain subgrid prescriptions. Unfortunately, yet unsurpris- 
ingly, a number of studies have clearly demonstrated that 
the results of these models are heavily affected by different 
choices for such prescriptions (e.g., Benson et al. 2003; Di 
Matteo et al. 2005), or for parameter values (e.g., Zavala et 
al. 2008). 

Thus extensive comparisons between different scenar- 
ios and data are generally conducted by means of semi- 
analytic modeling (SAMs) for baryons, often grafted onto 
gravity-only simulations for the dark matter (DM) evolu- 
tion. By definition of SAMs, the general behavior of the 
system is outlined a priori, and then translated into a set of 



(somewhat) physically-grounded analytical recipes — suit- 
able for numerical computation over cosmological timescales 
— for the processes which are thought to be more relevant 
to galaxy formation and evolution. 

Although SAMs should not be viewed as complete first- 
principles computations, they provide a convenient and pow- 
erful tool to test an assumed galaxy formation scenario (i.e., 
the general behavior and the adopted recipes) against exist- 
ing data, and to make predictions on future observations. 

In general, SAMs (e.g.. Cole et al. 2000; Hatton et al. 
2003; Cattaneo et al. 2005, 2006; Khochfar & Silk 2006; 
Baugh et al. 2005; Bower et al. 2006; Croton et al. 2006; 
Monaco et al. 2007; Somerville et al. 2008), apart from rel- 
atively minor variations, are constructed around two main 
assumptions: (i) the initial outcome of gas cooling within 
DM halos is, at any cosmic epoch, the development of a 
rotationally-supported disc (since Rees & Ostriker 1977; Silk 
1977; White & Rees 1978); these discs usually undergo mild 
to moderate star formation activity, unless extreme choices 
for the scaling of star formation efficiency with galaxy prop- 
erties are done; (ii) the most natural driver of episodes of 
violent star formation at any redshift is the merging of these 
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gas-rich discs, which in most models also constitutes the 
main channel for the formation of spheroids, and in partic- 
ular of large ellipticals (since Cole 1991). 

As a result of this disc-merger-driven framework, 
baryons tend to follow the hierarchical behavior of DM ha- 
los, and there is no inherent relationship between the mor- 
phology and the star formation history of galaxies. This is 
in sharp contrast with the basic observational fact that low- 
mass galaxies tend to be disc dominated, gas rich, blue, and 
actively star forming, whilst more massive galaxies tend to 
be red, gas poor, quiescent, and dominated by a spheroidal 
component mainly comprised of old stars. 

Due to these features, SAMs built around the two afore- 
mentioned assumptions — which from now on will be collec- 
tively referred to as 'standard SAMs' — tend to be in ten- 
sion with several observations (e.g., Somerville et al. 2008), 
manifested by the poor performances they had in anticipat- 
ing observational breakthroughs occurred more recently. For 
example, it is now well established that baryonic structures 
undergo the phenomenon referred to as 'cosmic downsizing', 
whereby massive star forming systems and associated super- 
massive BHs shined mostly at high redshift, while smaller 
objects display longer-lasting activity. Clearly, it is chal- 
lenging to obtain this behavior from the scheme outlined 
above; indeed, no model did until relatively recently, and 
the present situation remains unclear. In the past few years, 
almost all semianalytic teams introduced simple recipes of 
feedback from active galactic nuclei (AGNs) in their models, 
with the specific target of quenching star formation in high 
mass galaxies at low redshift. This additional ingredient sig- 
nificantly improves the situation, but does not directly alle- 
viate model difficulties in producing enough massive systems 
at high As a result, at least three state-of-the-art standard 
SAMs still do not correctly reproduce the downsizing trend 
in stellar mass, nor the archeological downsizing (Fontanot 
et al. 2009; see below for more details). 

A further example of challenges to models comes from 
the modest evolution of the cosmic star formation activity 
above z ~ 1 (the so called Madau plot), strikingly at vari- 
ance with model predictions (e.g.. Cole et al. 1994) generated 
before the advent of surveys effective in discovering dust- 
enshrouded star formation at high z (Madau et al. 1996, 
1998). It is fair to note that a fraction, but not all, of the 
discrepancy was due to the then adopted standard CDM cos- 
mology and a lower normalization for the fiuctuation spec- 
trum, resulting in significantly more rapid evolution at high 
redshift than the now favored ACDM. 

In addition, even latest and most refined SAMs are 
seriously challenged by the bright number counts and the 
high redshift peak of ^-distribution for submm galaxies. For 
instance, Baugh et al. (2005) showed that the only way 
to reproduce the statistic of submm sources, usually con- 
sidered the precursor of local ellipticals, in the context of 
their standard SAM, is to adopt an extremely top-heavy 
intial mass function (IMF) during galaxy-merger-induced 
starbursts. However, their model predicts masses of submm 
sources likely too low by more than one order of magnitude 
(Swinbank et al. 2008), and still shows discrepancies with 
observed trends of a/Fe in local ellipticals (Nagashima et 
al. 2005). 

Without doubt, the field of galaxy formation is led by 
observations. Indeed, physical processes have been continu- 



ously added to SAMs, or existing ones have been substan- 
tially revised by SAM developers in order to face serious mis- 
matches between model outputs and new datasets. Besides 
many relatively minor but subtle details, major examples 
comprise a treatment of the growth of supermassive BHs in 
galaxy centers and of the ensuing energetic feedback from 
nuclear activity (Granato et al. 2004 [G04j3; Bower et al. 
2006; Croton et al. 2006; Monaco et al. 2007; Somerville et 
al. 20080, the effects of 'cold' versus 'hot' accretion flows 
onto DM halos, as suggested by Dekel & Birnboim (2006) 
and implemented in a full SAM by Cattaneo et al. (2006; 
see also Somerville et al. 2008), or an extremely top-heavy 
flat IMF in merger-driven bursts (Baugh et al. 2005) . These 
examples show that the complexity and degrees of freedom 
of standard SAMs have been steadily increased by modelers 
in order to improve the agreement with the data, but de- 
spite these efforts several points of tension still remain (see 
Monaco et al. 2007). 

Within this paper, we follow a significantly different ap- 
proach and submit a novel scenario for galaxy formation, 
modifying the very basic assumptions of standard SAMs 
that, according to us, are the origin of their tensions with 
observations. Our scenario envisages that the fundamental 
dichotomy between the spheroid and disc components in 
galaxies reflects two dominating modes for the assembly of 
visible matter, feasibly being ultimately driven by two domi- 
nating modes governing the growth of DM halos (see below) . 
For typical L* galaxy halos, z ^ 2 corresponds to an era dom- 
inated by violent merging episodes, leading to huge bursts of 
star formation and to the observed coevolution of spheroids 
with hosted central supermassive BHs; 2^2 instead corre- 
spond to an era where the most relevant process is quiescent 
accretion of matter yielding, under suitable conditions, the 
formation of discs around preexisting spheroids. 

Many studies in the literature on the chemical and spec- 
trophotometric properties of local galaxy populations {stel- 
lar archeology; see Thomas et al. 2005 and references therein; 
Chiappini et al. 1997; Portinari & Chiosi 1999) reached the 
broad conclusion that galaxies of later type, which are less 
massive on average, formed their stars at later times and 
over a longer period (see also Gavazzi et al. 1996). This 
phenomenon is sometimes referred to as archeological down- 
sizing and is not reproduced by three state-of-the-art SAMs 
(Fontanot et al. 2009; however, see Kaviraj et al. 2005 for 
a discussion of color-magnitude relation in cluster ellipticals 
as a test for hierarchical models). These conclusions have 
more recently been conflrmed by modern surveys at high 
redshift, directly showing that the sites of active star forma- 
tion shift from high-mass galaxies at early times to lower- 
mass systems at later times [downsizing in time; Cowie et 
al. 1996; Guzman et al. 1997; Brinchmann & Ellis 2000; Ko- 
dama et al. 2004; Juneau et al. 2005; Bell et al. 2005; Noeske 
et al. 2007) . Further support for two different epochs and for- 
mation mechanisms of spheroids and discs comes from the 

^ in the context of a nonstandard SAM focused on the coevolu- 
tion of quasars and spheroids, of which this paper can be consid- 
ered an extension, see below. 

^ see also Hatton et al. (2003) and Cattaneo et al. (2007) for the 
effect of an highly idealized criterium of 'pseudo AGN feedback', 
phenomenologically inspired by the Magorrian (1998) relation- 
ship. 
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Figure 1. A schematic of tlie two-phase evolution of DM halos 
and associated modes of galaxy formation. 



analysis of the colour and structural properties of decom- 
posed galaxy components in the Millennium Galaxy Catalog 
(Driver et al. 2006). 

From the theoretical perspective, recent analysis of high 
resolution simulations of individual DM halos forming in cos- 
mological volumes (Zhao et al. 2003; Diemand et al. 2007; 
Hoffman et al. 2007; Ascasibar & Gottloeber 2008) have pro- 
vided support for our picture. In these studies, two distinct 
phases in the growth of DM halos have been identified: an 
early fast collapse featuring a few violent major mergers, 
and a later calmer phase including many minor mergers and 
smooth accretion. During the early collapse a substantial 
mass is gathered through major mergers, which effectively 
reconfigure the gravitational potential wells and cause the 
collisionless DM particles to undergo dynamical relaxation 
and isotropization (Lapi & Cavaliere 2009); therefrom the 
system emerges with a definite structure for the inner den- 
sity and gravitational potential (Lu et al. 2006) . During the 
later quiescent phase, moderate amounts of mass are slowly 
accreted mainly onto the halo outskirts, little affecting the 
inner structure and potential, but quiescently rescaling the 
overall mass upwards. Mo & Mao (2004) qualitatively sug- 
gested that this two-phase formation of DM halos may be 
at the origin of the disc-spheroid dichotomy, alleviating sev- 
eral problems of the standard SAMs. Here we take up this 
general idea, and construct a full semianalytic model capa- 
ble at making quantitative predictions to be compared with 
present and future observations. We explicitly point out that 
the backbone for the cosmological growth of DM halos is 
broadly the same as that adopted by all other SAMs. The 
novelty is in the fact that we identify the transition between 
the two phases, and we assume that the main processes driv- 
ing the evolution of luminous matter are strongly linked to 
the two different modes of DM assembly. 

Our model constitutes a natural extension to include 
disc formation at low-z, of the Antihierarchical Baryon Col- 
lapse scenario {ABC; G04; Granato et al. 2006; Silva et al. 
2005; Lapi et al. 2006; Mao et al. 2007; Lapi et al. 2008) 
that was focused on the high redshift coevolution between 
spheroidal galaxies and hosted supermassive BHs. This has 
been extremely successful in reproducing a wealth of obser- 



vations, including statistics of submm galaxies, properties of 
local elliptical galaxies, the results of deep K-hand surveys, 

nthe estimated local mass function of supermassive BH, and 
statistics of high-redshift quasars (QSOs). These results are 
essentially inherited by the model presented here. 

In summary, motivated by the successes of the high red- 
shift ABC framework and the compelling theoretical and ob- 

I servational evidence in support of a two-phase galaxy forma- 
tion scenario, we have developed a framework linking high 
and low redshift processes in order to generate the observed 
dichotomy between early-type and late-type galaxies. The 
plan of the paper is the following: we describe in detail the 
two-phase DM treatment in § 2; in § 3 we describe the mod- 
eling of the baryonic matter evolution for the spheroid and 
disc components; we present our results in § 4 and finally 
we summarize and discuss our findings in § 5. 

Throughout the paper we adopt the standard ACDM 
concordance cosmology, as constrained by WMAP 5-year 
data (Spergel et al. 2007). Specifically, we adopt a fiat cos- 
mology with density parameters flm = 0.27 and ~ 0.73, 
and a Hubble constant Ho — 70 km s~^ Mpc~^. 



2 DARK MATTER SECTOR 

In this work we compute the mass growth histories of DM 
halos using a binary mergertree with accretion based on the 
extended Press & Schechter formalism (see Lacey & Cole 
1993); specifically, we rely on the algorithm originally de- 
veloped by Cole et al. (2000), and recently improved by 
Parkinson et al. (2008) to reproduce the outcomes of A'^- 
body simulations. 

The algorithm starts from the expression for the mass 
fraction of a halo with mass M2 at redshift Z2 that was 
contained within a progenitor halo of mass Mi < M2 at 

Zl > Z2: 



/ = 



AS, 



(Acr2 



-(A5c) /2 Ao 



da^ 



AM 



(1) 



here AS, = Sc{zi)—5c{z2) is the difference between the linear 
thresholds for collapse at redshifts zi and Z2, while Aa^ — 
a^{Mi) — a^{M2) is the difference between the variances 
of the linear density fluctuations extrapolated at 2; = in 
spheres containing masses Mi and M2. 

Taking the limit of the above equation as zi —> Z2 one 
finds the merger rate as 



d/ 1 

d2 ^27r (Acr2)3 



dz 



da' 



dM 



(2) 



now one easily can work out the distribution for the number 
of halos with mass Mi into which a halo with mass AI2 splits 
during a step dz up in redshift: 



dN 



dMi dz Ml 



(3) 



Then given a mass resolution M^cs one may define the 
mean number of progenitors with masses Mi between Mres 
and M2/2: 



P = 



M2/2 



dN 
dAh 



dMi 



(4) 
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and the fraction of mass of the final object in progenitors 
below the resolution limit: 



diV Ml 



dMi M2 



(5) 



In fact, in constructing the above quantities we have re- 
placed 



dJV 

dMi 



AN 
dMi 



G{cri /(J2, 82 1 <J2) 



(6) 



where G{ai_/ 02, 82/(^2) is a perturbing function given by 
Parkinson et al. (2008), tuned to bring the Montecarlo 
merger histories in close agreement with the outcomes of 
state-of-the-art numerical simulations. 

Given all that, the rncrgcrtrcc algorithm is straightfor- 
ward: choose a mass M2 and redshift z for the final halo; 
pick up a step in redshift Az such that P 1 to ensure 
that the halo is unlikely to have more than two progenitors 
a.t z + dz; generate a uniform random number R between 
and 1; if J? > P do not split the main halo at this step, and 
simply reduce its mass to M2 {1 — F) to account for sub- 
resolution accretion; if 7? < P generate a random value of 
All between Afros and M2/2 consistent with the distribution 
given by Eq. (1), to produce two new halos with masses Mi 
and M2 (1 — F) — Mi; repeat the process on each new halo 
at successive redshift steps to build the overall merging tree. 

In our implementation of the algorithm we use an adap- 
tive step size Az such that P remains at a value significantly 
below unity, and then postprocess the tree by sampling it 
over convenient redshift intervals. In addition, we take the 
resolution mass Mics{z) as the one corresponding to a virial 
temperature of 10'* K, above which atomic gas cooling allows 
baryonic structures to condense. 

As an input of the algorithm, we use the Bardeen et al. 
(1986) power spectrum of density fluctuations with correc- 
tion for baryons (Sugiyama 1995), normalized so as to yield 
a mass variance ag « 0.8 on a scale of 8 h^^ Mpc. As an out- 
put, we obtain many realizations (up to several thousands 
within conceivable computational times) of the mergertree 
for a given present mass Mo at z = 0; each realization lists 
all the progenitors of Mo at different rcdshifts and describes 
how and when these merge together. We generate trees for 
masses Mo spanning the range from 10^ to 10^* Mq in loga- 
rithmic increments, and follow the related merging histories 
down to the resolution mass. 

We find that the halo growth along a given evolutionary 
track occurs in two distinct phcises: an early violent collapse 
where rapid growth is enforced by major mergers among 
several massive clumps; and a late period of gentle mass ad- 
dition through calmer accretion (see Fig. 1) extending down 
to the present time. A similar behavior has been pointed out 
in a number of recent numerical simulations (Weschler et al. 
2002; Zhao et al. 2003; Diemand, et al. 2007; Hoffman et al. 
2007; Ascasibar & Gottloeber 2008), and has been analyzed 
in semianalytic studies of Montecarlo merging trees to ex- 
plore the origin of the structural properties of DM halos (see 
Lu et al. 2006; Li et al. 2007; Lapi & Cavaliere 2009). 

In our view these different evolutionary phases of DM 
halo growth should significantly affect the main physical pro- 
cesses regulating the evolution of the baryonic matter within 
them; in particular, we envisage the violent eaxly collapse 
phase to be associated with the formation of the spheroid 



and hosted supermassive BH, while the gentle late phase to 
be favorable for the stable growth of galaxy disc around the 
preexisting spheroid-BH structure. 

As to the halo spatial structure we assume the standard 
NFW (Navarro, Prenk & White 1997) density profile 

p{r) ^ Avi, c^9(c) 

pc 3 (r/fvir) (1 -I- cr/rvir)^ 



(7) 



here pc is the critical density, g{c) = [log(l + c) — c/(l -|- 
c)]~^ is a weak function of the 'concentration' parameter 
c, and Avir ^ 187r^ + 82 [Q.Miz) - 1] - 39 [^Miz] - if is 
the non-linear collapse threshold in terms of the evolved 
matter density parameter Q,m(z) = (1 -I- 2:)^/[J1m (1 + 
z)'^ -\- JIa]. In fact, A^-body experiments (Taylor & Navarro 
2001; Zhao et al. 2003; Diemand et al. 2007; Hoffmann et 
al. 2007; Ascasibar & Gootloeber 2008) show that the NFW 
profile is established during the early fast collapse phase, 
with concentration parameter Ct — A for z>^zt\ in the slow 
accretion phase for z'^zt the DM halo potential well retains 
its shape and Vs stays put while the overall size rvir of the 
system increases, to the effect of rising the concentration 
parameter, see below. 

Equipped with these notions, we effectively trace the 
redshift evolution in the tree of a DM halo mass M^ with 
present value Mo as follows (see Fig. 1). First of all, we 
compute the concentration co of the mass Mo according to 
the prescription by Maccio et al. (2007): 



log Co = 1.071 -0.098 



log 



Mo 
Mo 



- 12 



(8) 



we stress that our computation neglects the scatter of cq at 
fixed mass, and any dependence of cq itself on the details 
of the merging history (see Weclisler et al. 2002; Zhao et al. 
2003). In addition, for the sake of simplicity we disregard 
the influence of the baryons on the halo structure in terms 
of adiabatic contractions or expansions (see Blumenthal et 
al. 1986; Gnedin et al. 2004). 

Then during the late, slow accretion phase we take M^ 
as the mass of the main progenitor, i.e., we follow only the 
main branch of the mergertree. This is achieved by starting 
at z = and working toward higher rcdshifts, taking the 
most massive halo at each splitting (merger) event. We ne- 
glect the baryonic processes (in particular star formation) 
occurring in the other branches of the tree; in other words, 
we make the approximation that when matter in the mi- 
nor branches of the tree joins the main progenitor it is still 
pristine, unevolved by baryonic processes, and thus mainly 
in the form of a gaseous medium with primordial metallic- 
ity. This is conceivable since: these halos have masses close 
to the cooling mass, and thus have just grown from below 
this threshold; we expect the star formation efficiency to be 
relatively low there, since supernova (SN) feedback would 
be most effective in ejecting gas from such shallow potential 
wells; such small halos are likely to contain small galax- 
ies that would not alter significantly the properties of the 
galaxy hosted by the main progenitor if they happened to 
merge with it. 

For a quantitative analysis, we compute and illustrate 
in Figs. 2 and 3 the properties of mergers undergone by DM 
halos along their growth history during the slow accretion 
phase. We find that DM halos grow predominantly through 
small accretion events and that the majority of our halos 
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Figure 2. Probability distribution for the number of mergers 
undergone by DM halos during the slow accretion phase (results 
over 1000 realizations); various linestyles refer to different ranges 
for the relative mass addition AM/M. In general, major mergers 
between halos are very rare events during the slow phase. 




Figure 3. Probability distribution for the relative mass addition 
to halos during the slow accretion phase (results over 1000 realiza- 
tions); various linestyles refer to the largest and the second largest 
value of AM/M in the tree during the slow accretion phase. In 
general the dominant mechanism for growing halos in the slow 
phase is the accretion of small lumps. 



(about 80%) do not undergo a substantial merger event; we 
recall that conventionally a major merger is defined as one 
in which the added mass exceeds that of the merging units 
by 1/3 or more, i.e. AM/M > 1/3, see Lacey & Cole (1993). 
Note, however, that a minority of halos (20% or so) do un- 
dergo a major merger event; plainly, within these systems 
the growth of stable galaxy discs can be temporarily inter- 
rupted. In addition, we note that the number of major merg- 
ers during the slow accretion phase is closely independent of 
halo mass; thus the average major merger rate is higher for 
more massive halos, which have more recently made their 
transition into the slow accretion phase. 

We compute the concentration Cz associated to the mass 
Mz after Zhao et al. (2003) using 



Figure 4. Probability distribution for the mass rank of the main 
progenitor in the mergertree at zt- In general, the main progenitor 
is actually the most massive halo at the transition. 



?.0 




0.0 0.5 1.0 1,5 2,0 

Mr.2/M„p 

Figure 5. Probability distribution for the mass ratio between the 
main progenitor and the most massive of all other halos in the 
mergertree at zt; in general, all the other halos are significantly 
less massive than the main progenitor. 



[log(l + e,)-c./(l + c.)]er^" = (!kV" (MlV'" ■ (9) 
[log(l + Co) - co/(l + Co)] Co \HoJ \MoJ 

here a — 0.48 (0.64) in the slow (fast) phase is a fitting 
parameter derived from TV-body simulations, and = 
Ho [^M (1 + z)^ + JIa]^^^ is the evolved 'Hubble constant'; 
a less accurate but simpler approximation of the expression 
above in the slow phase reads = co [Hq/ H^]^^^ with ~ 1, 
see Mo & Mao (2004). Following Zhao et al. (2003), the tran- 
sition between the slow accretion and the fast major merger 
phase occurs at the redshift zt where the concentration Cz 
decrease below the critical value Ct = 4. 

In Figs. 4 and 5 we illustrate the properties of the merg- 
ertree at Zt- First, we note that the main progenitor is gen- 
erally the most massive halo in the mergertree at zt ; second, 
the ratio between the mass of the main progenitor halo to 
the second most massive halo is generally very large. These 
findings show that, during the late slow accretion phase, the 
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Figure 6. Probability distribution for the transition redshifts 
(results over 1000 realizations); various linestyles are for differ- 
ent ranges of present halo mass Mq. The distributions of zt are 
roughly lognormal with a low redshift peak and a rather extended 
tail at high z. Note that lower mass halos have a higher zt on the 
average. 



Figure 8. Typical mass accretion history for a Milky Way-sized 
DM halo (Afo ~ 5 X 10^^ A/q). The transition between slow ac- 
cretion and fast collapse, that occurs at a redshift zt ~ l,is high- 
lighted as a dashed line. The halo growth (solid line) follows the 
main progenitor's one (dot-dashed line) during the late slow ac- 
cretion phase, and the composite halo (see text for details) during 
the early fast collapse. 
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Figure 7. Probability distribution for the ratio between the 
merging time and the dynamical time of the individual merg- 
ing units during the early, fast collapse phase (results over 1000 
realizations). In general, during the early collapse mergers occur 
on very short timescales. 



relevant evolution of the main progenitor is characterized by 
minor merger events with low mass halos. 

It is interesting to note (see Fig. 6) that on average 
the transition redshift zt decreases with increasing present 
day halo mass. In our scenario this means that larger halos 
have less or no time to develop a substantial disc compo- 
nent, hence massive halos will tend to host pure spheroids 
and viceversa, in broad agreement with observations. Inci- 
dentally, note that although massive halos spend the major- 
ity (or all) of their lifetimes within the fast collapse phase, 
the timescales for baryonic evolution, that are governed by 
the coevolution of the spheroid and the central supermassive 
BH, can be much shorter; in fact, these systems typically be- 
come 'red and dead' due to QSO feedback at relatively high 



redshifts, and passively evolve thereafter until the present 
(see § 3.1). 

At z > Zt, during the fast collapse phase, we compute 
Mz not as the main progenitor mass, but as that of the com- 
posite halo made of the overall mass in all the branches of 
the tree that will contribute to the main progenitor mass at 
Zt- Moreover, we consider only halos whose mass exceeds the 
critical halo mass for efficient gas cooling, i.e., we consider 
only halos where the related virial temperature is above lO"* 
K. 

Quantitatively, this can be justified by comparing the 
dynamical time for the composite halo with the major 
merger timescale. The former is defined as the timescale for 
the material to ballistically collapse to its center: 

in terms of the average density of the composite halo; the 
latter is defined as the average timescale for a major merger 
throughout the duration of the fast collapse phase: 

_ Nmm i^^s 

Tmorgc — r. ' 1-'--'-/ 

Olfast 

where Nmm is the number of times the main progenitor 
undergoes a major merger during the fast collapse phase, 
and 5ffast is the timelapse the halo spends within it. 

Fig. 7 shows that during the fast collapse phase several 
major merger events occur over timescales typically shorter 
than the dynamical time of the composite system, so just 
following the composite halo is a conceivable approximation. 

In Fig. 8 we plot one typical realization of the mass ac- 
cretion history for a Milky Way-sized DM halo, highlighting 
the transition redshift zt, the behavior of the composite halo 
at z > Zt and of the main progenitor at z < zt . 

Finally, we specify the angular momentum of the DM 
halo in the mergertree as follows. First of all, we recall that 
the angular momentum J is usually expressed in terms of 
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the dimensionless spin parameter A = J M~^^^, 
where E is the total energy of the halo; A'^-body experi- 
ments have shown that A does not correlate with halo mass 
nor concentration, is nearly independent of the redshift, and 
follows a lognormal distribution with average value 0.04 and 
scatter 0.5 dex (Bullock et al. 2001; Maccio et al. 2007). 

Thus for each halo we randomly select a value of A from 
such a distribution and neglect its evolution (e.g., Barnes & 
Efstathiou 1987; Kravstov et al. 1997; Vitvitska et al. 2002; 
Hetznecker & Burkert 2006). We are aware that the latter 
could impact on the evolution of disc properties, but choose 
to keep our treatment of this effect as simple as possible, 
given that to our knowledge a robust modeling has not yet 
been included into state-of-the-art SAMs; e.g., Somerville et 
al. (2008) use the spin parameter of the more massive halo 
at any given merger event. 



3 BARYONIC SECTOR 
3.1 Modeling the Spheroid 

The treatment of the baryonic processes in the fast major 
merger phase follows the recipes adopted by G04 to model 
the coevolution of spheroids and supermassive BHs, and al- 
ready exploited by our team in several previous papers (Silva 
et al. 2005; Cirasuolo et al. 2005; Lapi et al. 2006; Granato et 
al. 2006; Mao et al. 2007; Lapi et al. 2008). Here we provide a 
qualitative summary of the model focusing on its distinctive 
features, and defer the reader to G04 for all the details. 

We recall from § 2 that during the fast collapse phase, 
a rapid sequence of major mergers build up a DM halo of 
mass Mt at the transition redshift Zt; as for baryonic matter, 
we assume that condensation and cooling processes become 
effective at a formation redshift Zf > zt when the mass of 
the composite halo surpasses a substantial fraction of Aft, 
namely Mt/2 as widely adopted in the literature to define 
the 'formation epoch' (Lacey & Cole 1993; Kitayama & Suto 
1996). The results presented in this paper are almost insen- 
sitive to the exact choice for the fraction of Aft adopted 
in the definition of Zf. Specifically, we have checked that 
galaxy properties vary less than 10% when the mass frac- 
tion is changed between approximately 1/5 and 2/3, due to 
the strong effects of QSO feedback in terminating the star 
formation soon after Zf. 

After Zf, a mass A/inf « fbMt of baryonic matter, in 
cosmic proportion fb ~ 0.17 with the DM's, is shock heated 
to the virial temperature by falling into the gravitational 
potential well. This hot gas, assumed to follow an isother- 
mal distribution, may cool quickly especially in the denser 
central regions at the rate 



Afc, 



^cool 



in terms of the local cooling timescale 



^rool — 



3 pgas kT 



2fimpCni{r)A{T) 



(12) 



(13) 



here Pgas is the gas density, Ue is the electron density, T is 
the temperature, A(T) is the cooling function, and C ~ 5—10 
is a parameter describing the dumpiness of the gas. 

The cooled gas mass AfcoW, assumed to still follow the 



DM radial distribution, may form stars directly over the 
local dynamical timescale, providing a rate of star formation 



m = 

with 

^dyn 



dAfcoid(r,t) 
dr 



dr , 



37r 



32 Gp{r) 



1/2 



(14) 



(15) 



This is conceivable during this evolutionary stage since the 
ongoing major mergers continuously reshuffle the gravita- 
tional potential, enforcing dynamical relaxation and orbit 
isotropization of the coUisionless DM and stellar components 
(Lapi & Cavaliere 2009). We recall that usually SAMs as- 
sume instead that the first result of gas cooling is the forma- 
tion of rotationally supported discs, characterized by much 
milder star formation activity, since the adopted star forma- 
tion timescale is typically much longer than some dynamical 
timefl 

With our prescriptions, large galactic halos can attain 
star formation rates of the order of ~ 1000 solar masses per 
year over timescales of a fraction of Gyr. This is required 
to explain the submm galaxy population without invoking 
an extremely top-heavy IMF (e.g., Baugh et al. 2005). In 
fact, our IMF has the standard Salpeter slope 1.25 in the 
high mass tail, and flattens to a slope 0.4 below IMq. As 
shown in Romano et al. (2005), this performs better than the 
Salpeter one in reproducing the detailed chemical properties 
of elliptical galaxies. 

Star formation promotes the gathering of some cool gas 
into a low-angular-momentum reservoir around the central 
supermassive BH. A viable mechanism for this process is ra- 
diation drag (see discussion by Umemura 2001; Kawakatu & 
Umemura 2002; Kawakatu, Umemura & Mori 2003), which 
has the nice feature of predicting a mass transfer rate to the 
reservoir proportional to the SFR to a good approximation: 



A/rd = Qrd X 10 



(16) 



The constant of proportionality ord ~ 1 — 5 can be fixed 
to produce a good match to the correlation between the 
spheroid and the supermassive BH masses observed in the 
local universe. The quantity 



TRD 



TO 



z 



1012 M, 



1/3 



(17) 



represents the effective optical depth of the gas clouds in 
terms of the normalization parameter tq ~ 1 — 5; for more 
details, see the discussion around Eqs. (14) to (17) in G04. 

Eventually, this gas accretes onto the BH powering the 
nuclear activity; in this early phase plenty of material is 
supplied to the BH, so that the latter can accrete close to 
the Eddington limit 



A^BH = AEdd 



1 - V A/bh 



(18) 



and grows almost exponentially from a seed of 10 Mq . The 
e-folding time involves the Eddington time t-Edd ~ 4 x 10* yr, 

for instance, SOtdyn in Hatton et al. (2003), ~ 200tdyn in Cole 
et al. (2000), ~ 350tdyn in Bower et al. (2006), and ~ 15tdyn in 
Croton et al. (2006). 
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the radiative efficiency 77 ~ 0.15, and the actual Eddington 
ratio AEdd ~ 0.3 — 3. 

The energy fed back to the gas by SN explosions and 
BH activity regulates the ongoing star formation and BH 
growth. The two feedback processes have very different de- 
pendencies on halo mass and on galaxy age (e.g., on the time 
since Zf). The feedback due to SN explosions removes the 
starforming gas at a rate 



Msn = -- esN — i]{t) 



(19) 



here is the velocity dispersion within the bulge, iSsN ~ 
10*^ erg is the energy released in a single SN event, 77SN 
is the number of Type II SNe expected per solar mass of 
formed stars (determined by the IMF), and esN ~ 0.05 is 
the fraction of this energy which is effectively coupled to the 
gas. Thus the SN feedback evolves almost in parallel with the 
star formation; it is very effective in low-mass halos severely 
limiting the growth of stellar and BH components there, but 
is of minor importance in the more massive galactic halos. 

The QSO feedback considered by G04 acts both on the 
cold as well as on the hot gas, unbinding them from the DM 
halo potential well at a rate 



Mqso 



-2 X 10'' 



L 



3/2 



eqso 



(cr/300kms- 



(20) 



this functional form is suggested by theoretical models of 
line-driven winds and observations of BAL QSOs (see the 
derivation leading to Eqs. [29] to [31] in G04). The Edding- 
ton luminosity I/Edd,46 5 in units of 10*® erg s~^ is a conve- 
nient measure of the BH mass, and eqso ~ 1 — 5 is a strength 
parameter. 

As a consequence, the QSO feedback grows exponen- 
tially during the early phases of galaxy evolution, following 
the exponential growth of the supermassive BH mass. It is 
is negligible in the first 0.5 Gyr in all halos, but abruptly 
becomes notably important in DM halos more massive than 
10^^ Mq^ structures weakly affected by SN feedback. Even- 
tually, in these systems most of the gas becomes unbound 
from the potential well of the galaxy halo (see Lapi, Cava- 
liere & Menci 2005 for the impact of QSO feedback on galaxy 
groups and clusters), so that star formation and BH activ- 
ity itself comes to an end on a timescale which is shorter for 
more massive galaxies. 

Indeed, the positive feedback on BH growth caused by 
star formation, in cooperation with the immediate and neg- 
ative feedback of SN, and the abrupt and dramatic effect of 
QSO feedback, are able to reverse the formation sequence 
of the baryonic component of galaxies compared to that of 
DM halos; the star formation and the buildup of central BHs 
are completed more rapidly in the more massive halos, thus 
accounting for the phenomenon now commonly referred to 
as downsizing. 

Before QSO feedback dominates the evolution, radia- 
tion is highly obscured by the surrounding dust. In fact, 
these protogalaxies are extremely faint in the UV-optical 
rest frame and are more easily selected at submm wave- 
lengths. The nuclear emission is also heavily obscured, and 
easier to detect in the hard X-ray band. On the other hand, 
when the central supermassive BH is massive and powerful 
enough to remove most of the gas and dust from the sur- 
roundings, the active nucleus shines as an optical QSO. Fol- 
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Figure 9. A schematic of the disc formation geometry. 

lowing this stage, the BH is already present at the galaxy 
center, thus any subsequent supply of gas to the spheroid 
produces an immediate QSO feedback, and thus is unable 
to substantially affect the stellar or BH mass: afterwards, 
the stellar populations in the spheroid evolve largely in a 
passive manner. 

Other SAMs (e.g.. Bower et al. 2006; Croton et al. 2006) 
introduced the 'radiomode feedback', which is active only 
in massive objects and at late times to halt cooling flows, 
but has no effect during the principal growth phase of most 
galaxies and AGNs. Also the highly idealized pseudo AGN 
feedback considered by Hatton et al. (2003) and Cattaneo et 
al. (2007) Q is somewhat representative of radiomode feed- 
back. By converse, in G04 a central role is given to the possi- 
ble feedback originated by the main episode of supermassive 
BH growth, which is responsible for the QSO activity at high 
z. 

The model described above has proved to be extremely 
successful in reproducing a wealth of observations, includ- 
ing statistics of submm galaxies, properties of local ellip- 
ticals, the results of deep Tf-band surveys, demography of 
supermassive BH relics, and statistics of high redshift QSOs. 
These successes are essentially inherited by its generalization 
presented here; in fact, we keep the model parameters fixed 
to the values used in the papers by Lapi et al. (2006) and 
Mao et al. (2007). We list the model parameters and their 
fiducial values in Table 1, stressing their relative relevance 
in the present context. 



3.2 Modeling the Disc 

At z < zt during the slow accretion phase, conditions 
become sufficiently quiescent to allow the dissipationless 



* They simply stops cooling when ^ Af^uig 
the sum is over all the galaxies in a halo. 
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Table 1. Model Parameters 



Description Symbol Fiducial value Reference in the text Impact on this work 



Spheroid + BH (ABC) 



Clu-inping fcLctor 


c 


7 






R,£Ldi3.tion drcig efficiency 




9 ^ 




Mild 


TNI ("irvci n 1 1 '7Ci't"i nvi rif rrntir'til rif^T^tn 


To 


I 




Weak 


BH radiative efficiency 


V 


0.15 


Eq. (18) 


Mild 


Eddington ratio 




1 


Eq. (18) 


Weak 


SN feedback efficiency 


eSN 


0.05 


Eq. (19) 


Strong 


QSO feedback efficiency 


fQSO 


1.3 


Eq. (20) 


Strong 


Disc (vdBOl) 










Star formation efficiency 


fsf 


2.5 X 10"* 


Eq. (24) 


Strong 


Schmidt law exponent 


n 


1.4 


Eq. (24) 


Mild 


Gas velocity dispersion in Toomre 




6 km s" 


Eq. (26) 


Weak 


Normalization constant in Toomre 


Q 


1.5 


Eq. (26) 


Weak 


SN feedback efficiency 




10-4 


Eq. (27) 


Strong 


Dust (GRASIL) 










Fraction of gas in molecular clouds 


/mc 


0.25 


Sect. 4.3 


Mild 


Optical depth of molecular clouds (at 1 /im) 


tmc 


60 


Sect. 4.3 


Weak 


Escape time from molecular clouds 


tesc 


2.0 Gyr 


Sect. 4.3 


Weak 



Note. - A Romano IMF </>(m*) is adopted: (pim*) oc m^, ^'^^ for m* > Mq and (f>{mi,) oc rra^ '^ ^ for m* < Mq. 



growth of discs from accreting materiaJ. Depending on zt 
and the shape of the individual growth history, under suit- 
able circumstances a substantial disc component may de- 
velop. To describe the process, we adopt the model by van 
den Bosch (2001 [vdBOl]); here we provide a quick overview 
of it, but defer the reader to the original paper for its full 
description. 

At each time step new baryons are accreted onto the 
halo at a rate Mz proportional to the DM's, in terms 
of the universal baryon to DM fraction « 0.17. As this 
material enters the halo it is assumed to be heated to the 
virial temperature, and to be distributed with an isothermal 
profile. The angular momentum distribution of the hot gas 
mirrors that of the DM component, so that the change in 
the angular momentum over a time interval St reads 



5 J = J{t) - J{t - 5t) , 



(21) 



where the halo total angular momentum J = 
G M^^^ X/\E\^^'^ is specified in terms of the halo spin 
parameter A, see § 2. 

Equating the gained angular momentum to that of a 
uniformly rotating shell of material, one obtains the circular 
frequency (see also Fig. 9): 



UJO 



r p{r) dr 



(22) 



The gas is then allowed to cool and collapse, conserving the 
initial angular momentum gained from the DM halo. The 
timescale for condensation icoii = max(tdyn, icooi) is given 
by the maximum of the dynamical and cooling time. After 
a time t' = t + tcow, the cooled gas is added to the disc 
annuli with radius corresponding to where it becomes 
centrifugally supported upon dissipationless collapse from 
the original cylindrical shell radius Ri, i.e.. 



Ri = 



Wo 



1/2 



(23) 



in terms of the local circular velocity 14 (r) = [GM(< 
r)/r]^^^. Thus the disc is allowed to grow in an onion- like 
fashion, and in this computation no specific disc profile is 
adopted a priori. 

The star formation rate is then assumed to follow the 
empirical Schmidt (1959) law, i.e., it is related to the surface 
density E(r, t) of cold gas in the disc through: 



V'(r,t) = esf 



E(r,i) 



Mq kpc 



Mq kpc yr 



(24) 



here e^f ~ 2.5 x 10^^ is a fudge parameter controlling the 
star formation efHciency, and n ~ 1.4 is fixed to match the 
properties of local spiral galaxies (Kennicutt 1998). 

At each time step and for each annulus within the disc, 
we compute the amount of material converted from gas to 
stars by solving 

dE(r,t) 



dt 



= -V(r,t) 



(25) 



actually we also impose that a gaseous disc annuli becomes 

eligible for star formation only once its surface densities sur- 
passes a critical threshold Scrit given by the Toomre (1964) 
criterion: 

where k{R) is the epicycle frequency (see vdBOl for details), 
Q ~ 1.5 is a constant, and cTgas ~ 6 km s"^ is the velocity 
dispersion of the gas. 

Cool gas may be removed from the disc through SN 
winds; we compute the related mass depletion in a way anal- 
ogous to the spheroidal modeling: 



MsN(r,i) 



2£SN ?7SN -EsN 



'4>{r,t) ■ 



(27) 



Vclc(r,t) 

here in the denominator the local escape velocity Vesc is used. 
Finally, we model the chemical evolution of the stellar mar 
terial on using the instantaneous recycling approximation. 
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Figure 10. Redshift evolution of the baryonic mass components 
for a Milky Way sized galaxy at z = 0: DM mass (solid/black), 
disc gas (dotted/blue), disc stars (short-dashcd/green), bulge gas 
(dot-dashed/red), bulge stars (long-dashed/orangc). The two ver- 
tical lines mark the transition and formation redshifts. 



Figure 11. Decomposition of the overall rotation curve (solid 
line) for our Milky Way-type galaxy at z = in terms of the con- 
tributions from DM (dashed/green line), bulge (dot-dashed/blue 
line) and disc (dotted/blue line). 



For the sake of coherence, and at variance with vdBOl, 
in the disc modeling we adopt the same IMF used in the 
treatment of the spheroid evolution. We find that the re- 
sults concerning the disc structure are affected by 10% from 
other reasonable choices of the IMF; we may recover reliable 
matches to the properties of the local galaxy population with 
all commonly used IMF by altering the fudge parameters of 
Table 1 within their physical limits. 

The main differences in our modeling with respect to 
vdBOl are the following: we use individual growth history 
derived from our detailed mergertree, while vdBOl adopts 
only an averaged smooth fit; we use a the prescription by 
Maccio et al. (2007) for the concentration parameter of DM 
halos, while vdBOl relies on Bullock et al. (2001); finally, at 
variance with vdBOl we take into account the gravitational 
effect of the preexisting spheroid on the dynamics of the 
forming disc but neglect to model the adiabatic response of 
the DM halos to disc settling. 



4 RESULTS 

4.1 Disc structure and dynamics 

In this Section we analyze the behavior of individual galax- 
ies, focusing on the disc component; the formation and evo- 
lution of the spheroidal component has been extensively con- 
sidered in several previous papers by our team (G04, Silva 
et al. 2005, Cirasuolo et al. 2005, Lapi et al. 2006). 

To this purpose we select a fiducial model galaxy with 
current mass Mo « 5 x 10^^ Mq similar to that of the 
Milky Way (see Xue et al. 2008; Naab & Ostriker 2006), 
and with an average spin parameter A = 0.04 (see Maccio 
et al. 2008). We find that the resulting galaxy components 
at z = are in generally good agreement with observed 
Milky Way properties, finding that Mdisc.o ~ 1-5 x Mq, 
Mbuige.o ~ 1.3 X 10^ Mq, see also Naab & Ostriker (2006). 

In Fig. 10 one can see the buildup of the various com- 
ponents for this fiducial galaxy. At 2: > zt a strong growth of 
the spheroidal component takes place which is halted by the 
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Figure 12. Decomposition of the disc surface density profiles 
(solid line) for a Milky Way sized galaxy at z = 0, in terms of the 
stellar (dotted-line) and gaseous (dashed line) components. 

QSO activity after approximately 10* years from z/. Fol- 
lowing this, the stellar populations in the spheroid evolve 
passively, and the residual gaseous material is originated 
from the stellar recycling. At z < zt new gas quiescently 
accretes onto the halo forming a disc structure; note that 
star formation is delayed until the cold gas surface density 
becomes sufficiently large to overcome the critical star for- 
mation threshold. 

In Fig. 11 we present the rotation curve decomposition 
for our fiducial galaxy. The total rotation curve is flat out 
to large radii; there it is DM dominated, while in the in- 
ner regions the baryonic components of the disc and the 
bulge dominate the gravitational potential, in agreement 
with kinematic models of the Milky Way; we flnd a peak 
rotation velocity VJnax ~ 210 km s~^. 

Fig. 12 shows the disc radial surface density profile at 
z = 0; we obtain general exponential stellar profile out to 
a truncation radius of 16.2 kpc, in broad agreement with 
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Figure 13. Redshift evolution of the overall surface density pro- 
file for a galaxy with spin parameter A = 0.06 resulting in an 
extension that is twice the fiducial model to highlight the evo- 
lution. Note that at 2 0.8 a minor merger event occurs, the 
exponential disc structure is temporarily disrupted, and an anti- 
truncated disc (see text for details) develops. 
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Figure 14. Impact of preexisting bulge masses on the disc surface 
density profiles (solid line) for a Milky Way-sized galaxy halo 
Mo ~ 5 X 10^^ Mq at 2: = 0; as the bulge mass approaches the 
disc mass of about 10^" Mq discs become more compact. Higher 
bulge masses yield significantly more concentrated discs. 

the value of 12 kpc observed for the Milky Way (see Naab 
&i Ostriker 2006). The gaseous disc is more extended than 
the stellar one due to the critical star formation thresh- 
old. However, note that the gaseous disc is depleted in the 
central regions because of star formation, and there the 
stellar component dominates. Fitting an exponential pro- 
file Eoe"'''''''' to the stellar disc, we find an exponential 
scale radius ra ~ 2.6 kpc, which is comparable to observa- 
tional estimates 2.5 — 3.5 kpc for the Milky Way, see Sackett 
(1997). Finally, within this model, we find a BH mass at 
2 = of Mbh,o ~ 6.7 X 10^ A/0, which is less massive but 
still consistent with the one at the center of the Milky Way 
Mbh.o ~ 3.6 X 10" Mq, see Eisenhauer et al. (2005). 

For the sake of completeness we highlight the related 
redshift evolution of the disc profile in Fig. 13. For illustra- 



tive purposes it is clearer to have an elongated disc profile; 
thus we choose a realization with a spin parameter A — 0.06 
larger than the fiducial value. The disc (gas and stars) natu- 
rally evolves from the inside-out, retaining a quasi exponen- 
tial profile with scalelength increasing over time. Although 
the vast majority of model discs have a quasi-exponential 
surface density profiles, the detailed shape depends on the 
details of the specific growth history, and in particular on 
the transition redshift zt. 

We stress that our discs develop following the buildup 
of the spheroidal inner component, which affects the over- 
all gravitational potential. Fig. 14 illustrates how the z — 
disc surface density depends on the mass of a preexisting 
bulge. We find that for identical realizations for our fidu- 
cial model, but imparting a bulge mass by hand, the disc 
structure (and thus its evolution) is altered. This dynamical 
interdependence results in a disc structure which becomes 
significantly more compact as the bulge mass approaches the 
disc mass at = 0. In addition, we find that with the pres- 
ence of a substantial spheroid component, higher transition 
redshifts zt yield more extended discs at z = 0. 

Recent observational studies have shown that in general 
exponential discs come in three categories, corresponding to 
simple exponential (type I), truncated (type II) and anti- 
truncated (type III) surface density profiles (Pohlen & Tru- 
jillo 2006). We find that exponential surface density profiles 
(of type I and II) are generated as a generic feature of our 
model (see Fig. 12). Interestingly, discs which have under- 
gone a recent minor merger event have a significant amount 
of material added to the outer parts in a non smooth fashion, 
and these mimic disc anti-truncations (type III, see Fig. 13 
at 2 ~ 0.8). However, the disc anti-truncations are typi- 
cally masked through subsequent accretion, that restore the 
overall exponential surface densities. We plan to address this 
issue in a subsequent work. 



4.2 Galaxy properties at z — 

Now we turn to study the properties of the local galaxy pop- 
ulation; to this purpose we generate catalogues of galaxies 
that encompass a representative range of 2 = halo masses, 
from 10^ to 10^** Mq in logarithmic increments. We then 
exploit the statistics of halos containing one single galaxy, 
namely, the galaxy halo mass function (GHMF) as provided 
by Shankar et al. (2006); the latter authors provide the fol- 
lowing analytic fit 

e(A/o) = -|(^)%-^^°/*, (28) 

with a ^ 1.84, 61 3.1 x 10"* Mpc-^ and M ^ 1.12 x 
10^'^ Mq . This function is derived by subtracting the group 
and cluster mass function (Martinez et al. 2002) from the 
Sheth & Tormen (1999) mass function. Thus, for halos with 
Mo < a few 10^^ Mq the GHMF closely follow the Sheth 
& Tormen mass function, while the falloff at larger masses 
mirrors the increasing probability of multiple occupation. In 
principle, we may account for the halo occupation number 
through the full mergertree, but this will introduce uncer- 
tainties related to poor knowledge of processes like dynami- 
cal friction, tidal stripping etc. etc.; thus we prefer to bypass 
this problem using the GHMF. 
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Figure 15. The relation between the host DM halo mass and the 
stellar mass ratio. Model results (over 1000 realizations) for indi- 
vidual galaxies (crosses), model average (solid line) and quartilos 
(dashed lines) are compared to the observational determination 
by Shankar et al. (2006, shaded area). 



Figure 17. Gas to stellar mass fraction as a function of the 
stellar mass. Model results (over 1000 realizations) for individ- 
ual galaxies (crosses), model average (solid line) and quartiles 
(dashed lines) are compared to the observational determination 
by Baldry et al. (2008, shaded area). 
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Figure 16. Fraction of baryons converted into stars as a function 
of host halo mass. Model results (over 1000 realizations) for indi- 
vidual galaxies (crosses), model average (solid line) and quartiles 
(dashed lines) are compared to the observational determination 
by Shankar et al. (2006, shaded area). 



Figure 18. The occurrence of bulge to total mass ratio in our 
model, binned for different halo mass ranges (results over 1000 
realizations); disc-dominated galaxies occur preferentially in low- 
mass halos, while spheroid-dominated galaxies occur preferen- 
tially in massive halos. 



In Figs. 14-16 we compare our model predictions with 
the results of Shankar et al. (2006) and Baldry et al. (2008), 
who derive a number of galaxy properties as a function of 
the host DM halo mass. In Fig. 15 we consider the frac- 
tion of stellar to total mass within DM halos; we see that 
there is a steep increase in the DM dominance for low mass 
halos, since these provide inefficient environments for star 
formation mostly due to the impact of SN feedback. 

Fig. 16, constituting a different rendition of the previ- 
ous plot, directly highlights the fraction of available baryons 
condensed into stars as a function of the host halo mass. 
In halos of masses exceeding few lO^'^ Mq star formation 
is more efficient; in the absence of a substantial impact of 
QSO feedback the efficiency would keep growing with in- 
creasing mass (despite an increasing difficulty of the cooling 



processes), while both the data and our model show a clear 
flattening. 

Fig. 17 illustrates the correlation between gas fraction 
(ratio of the total cold gas to the total baryonic mass within 
the galaxy) and the overall stellar mass, compared to the 
data by Baldry et al. (2008); less massive galaxies typically 
have a significantly larger gas fraction. This is because more 
massive galaxies are typically spheroid dominated and thus 
underwent strong gas ejection by the QSO feedback. On 
the other hand, lower mass galaxies are typically disc dom- 
inated, QSO feedback thus is relatively unimportant, and 
also the critical surface density threshold becomes increas- 
ingly difficult to surpass. 

Fig. 18 illustrates the occurrence of bulge to total mass 
ratio in our model, binned in host halo mass. The behavior 
of the galaxies in our model is dichotomic, with the disc- 
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log(M„,o) [Mq] 

Figure 19. _B-band luminosity as a function of the host DM 
halo mass. Model results (over 1000 realizations) for individual 
galaxies (crosses), model average (solid line) and quartiles (dashed 
lines) are compared to the observational determinations collected 
by Tonini et al. (2006, shaded area). 
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Figure 20. The TuUy-Fisher relation. Model results (over 1000 
realizations) for individual galaxies (crosses) , model average (solid 
line) and quartiles (dashed lines) are compared to the observa- 
tional determination by Giovanelli et al. (1997, shaded area). The 
Possinian error bars illustrate the relative abundance of galaxies 
with different circular velocity. 

dominated galaxies to occur preferentially in low-mass hales, 
while spheroid-dominated galaxies to occur preferentially in 
massive halos. This result is basically linked to the distribu- 
tion of transition redshift zt discussed in § 2 and illustrated 
in Fig. 6. 

Although not reported here, we stress again that this 
model inherits from the ABC scenario the good match with 
the observed local BH mass vs. bulge relationships (see G04, 
Cirasuolo et al. 2005, Lapi et al. 2006). 

4.3 Spectrophotometric properties 

In order to analyze the luminous properties of galaxies, 
we interface our model with the spectrophotometric code 



GRASIL (Silva et al. 1998), that accounts for the attenua- 
tion and reradiation of starlight by dust. 

GRASIL uses stellar population synthesis models based 
on the Padova evolutionary tracks, which include the effects 
of dusty envelopes around asymptotic giant branch stars 
(Bressan, Granato & Silva 1998). Then each single stellar 
population is summed taking into account the appropriate 
age and metallicity, and weighted with the star formation 
rate to obtain the unattenuated SED 

Fx{r)= [ (;xlr~t,Z{t)]mdt; (29) 
Jo 

here r is the age of the galaxy, t is the birth age of an 
individual single stellar population (\{t, Z), and xp{t) is the 
star formation rate. 

For the detailed description of dust attenuation and re- 
processing of starlight we defer the reader to the papers by 
Silva et al. (1998, 2005). We just stress here that GRASIL 
includes the effect of differential dust extinction of stellar 
population, i.e., younger stellar generation are more affected 
by dust obscuration; this is because stars form in molecular 
clouds, an environment denser than the average, and pro- 
gressively get rid of them. 

The GRASIL SEDs depend on the following basic pa- 
rameters: the fraction /mc of gas in the form of molecular 
clouds rather than in the diffuse interstellar medium; the op- 
tical depth tmc of molecular clouds to the radiation emitted 
from a source at their center (at 1 fim); the escape time tasc 
of newly born stars from molecular clouds. On the basis of 
previous works in which GRASIL has been coupled with var- 
ious SAMs, and in particular with the one of the Durham 
team (see Granato et al. 2000 for details; also Baugh et al. 
2005, Monaco et al. 2007), we set the GRASIL parameters 
to the standard values reported by Silva et al. (2005) and 
listed in Table 1. Note that since in this paper we do not 
consider regions of the SED strongly affected by dust emis- 
sion, the dependence of our results on these parameters is 
mild/weak. Due to the preliminary nature of this work, we 
do not exploit the multiwavelength capabilities of GRASIL 
to the full extent, but focus on reproducing several local 
galaxy population properties in selected bands; we delay a 
more refined analysis for future work. 

In Fig. 19 we show the _B-band luminosity of our model 
galaxies as a function of DM halo mass. We find a relatively 
strong correlation with little scatter, in general agreement 
with the data by Tonini et al. (2006). The break around 
10^^ Mq is due to the impact of SN feedback in small sys- 
tems, where star formation becomes progressively less effi- 
cient. 

Fig. 20 illustrates the I-band TuUy-Fisher relation from 
our model; note that we extract the peak rotation velocity 
from our model by fitting an exponential to the disc sur- 
face density profile. Then we use the obtained scale radius 
rd to define the maximum velocity as Vmax ~ 14(2.2 r^) 
in analogy with the observational methods; this procedure 
therefore does not resort to further assumptions about disc 
structure and dynamics. 

Our model result is compared with the data by Gio- 
vanelli et al. (1997), finding excellent agreement in both 
slope and normalization; fitting our result with the law 
Mi = m [logj^Q(V2.2) — 2.5] + c, we obtain m — —8.30, and 
c = —21.19, to be compared with the observational values 
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m — —7.68 and c = —21.0. We also represent the intrinsic 
scatter in our results by the blue dashed contours, showing 
a slight if systematic increase in scatter towards higher ro- 
tational velocities, in general agreement with observations 
(Giovanelli et al. 1997). Finally, we show the relative abun- 
dances of galaxies as represented by the Poissonian error 
bars that account for the relative numbers of halos within a 
cosmological volume. We see that the slowly rotating, faint 
galaxies are more abundant; galaxies become rarer as we 
move toward larger circular velocities, and within our sam- 
ple we do not find any galaxies with rotational velocities 
larger than V2.2 > lO^'® km s"\ 

In Fig. 21 we present the r*-band luminosity function 
from our model at z = 0, and compare it to the fit by Benson 
et al. (2007) based on SDSS data; we find an overall good 
agreement. We confirm, in tune with a number of previous 
works, that the fiattening at the faint end is mainly due to 
the impact of SN feedback in small systems (Benson et al. 
2003 and references therein), and that the steepening at the 
bright end is mainly caused by the impact of QSO feedback 
in massive galaxies (Somerville et al. 2008 and references 
therein). In addition, we highlight the different contribu- 
tions to the overall luminosity function from the spheroid 
and the disc component; the latter typically dominate the 
faint end, while the former dominate the bright end, as ex- 
pected on an intuitive basis (see also Kauffmann, White & 
Guiderdoni 1993; Baugh, Cole & Frenk 1996; Kauffmann 
& Chariot 1998; Somerville & Primack 1999; Hatton et al. 
2003; Tasca & White 2005). 

We stress that the simultaneous fitting of the galaxy 
luminosity function and of Tully-Fisher relation is challeng- 
ing for many SAMs, and so constitutes a big success of our 
model (see Courteau et al. 2007; Dutton et al. 2007; Bell et 
al. 2003). 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have proposed a novel scenario for the for- 
mation and evolution of galaxies in the standard ACDM 
framework. 

We have been motivated by several recent high- 
resolution A'^-body simulations (Zhao et al. 2003; Diemand 
et al. 2007; Hoffmann et al. 2007; Ascasibar & Gottloeber 
2008), that recognize the DM halo growth to occur in two 
rather distinct phases: an early violent collapse featuring 
a few major mergers, and a late quiescent accretion onto 
the halo outskirts that does not affect the inner regions 
where the galactic structure resides. We associate these two 
phases to two different modes of galaxy formation, leading 
to spheroids and discs. 

Specifically, we envisage that spheroids form during the 
fast collapse phase, when violent major mergers reshuffle 
the gravitational potential and cause dynamical relaxation 
and orbit isotropization of the DM and stellar components 
(see Lapi & Cavaliere 2009). Meanwhile, strong starburst 
activity and the growth of a central supermassive BH take 
place in parallel. The ensuing SN explosions and the nuclear 
activity feed energy back to the baryons, and regulate the 
ongoing star formation rate and BH growth. These mutual 
energy feedbacks actually reverse the formation sequence of 
the baryonic component of galaxies compared to that of DM 
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Figure 21. The r*— band luminosity function. The overall results 
(over 1000 realizations), with the contributions from spheroids 
(red hatched region - upper panel) and discs (blue hatched region 
- lower panel) highlighted, is compared to the fits derived from 
the SDSS data by Benson et al. (2007, coloured dashed lines). 



haloes: the star formation and the buildup of central BHs 
are completed more rapidly in the more massive haloes, thus 
accounting for the phenomenon now commonly referred to as 
downsizing. In the subsequent slow accretion phase, during 
which major mergers are rare, the quiescent growth of a 
disclike structure around the preformed spheroids can occur 
by dissipantionless collapse. 

We then test this new scenario against observations 
resorting to the semianalytic technique. To this purpose, 
we adopt standard and widespread models. As to the DM 
evolution we base on the algorithm by Cole et al. (2000) 
and Parkinson et al. (2008) supplemented by the results of 
the A-body experiments by Zhao et al. (2003). As to the 
spheroid component we rely on the prescriptions by G04 
and following developments. As to the disc buildup we base 
on the recipes by vdBOl. Finally, we couple everything to 
the spectrophotometric code GRASIL by Silva et al. (1998). 

Note that the SAM developed here can be viewed 
as an extension of our previously proposed ABC model 
(see G04). The latter dealt with the high-redshift spheroid- 
supermassive BH formation, and proven to be successfuU in 
many respects (see Silva et al. 2005, Cirasuolo et al. 2005, 
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Lapi et al. 2006, Mao et al. 2007, Lapi et al. 2008). Prac- 
tically, we now include the disc formation at low redshift, 
so extending it to encompass all morphological galaxy types 
and cosmic epochs. 

Though we are confident to have described through con- 
ceivable physical recipes the key processes ruling galaxy for- 
mation and evolution, we must admit that our modeling 
disregard or treat crudely several aspects that may play a 
relevant role: baryon impacts on the detailed structure of 
DM halos, environmental effects, angular momentum evo- 
lution, bar instabilities, halo occupation distributions, etc. 
In fact, in this preliminary study we pursue the strategy of 
'keeping it as simple as possible', our aim being to test with 
minimal ingredients whether our scenario could provide re- 
sults in accord at least with the local galactic observables, 
and eventually it proved to perform such a remarkable task 
surprisingly well. 

Specifically, we have shown our model to reproduce the 
observed stellar mass fractions (see Figs. 15 and 16), gas 
content (see Fig. 17), morphological dichotomy (see Fig. 18), 
mass-to- light ratios (see Fig. 19), TuUy-Fisher relation (see 
Fig. 20), and luminosity functions (see Fig. 21) of the lo- 
cal galaxy populations. In future works we aim to compare 
our model predictions to the intermediate- and high-redshift 
data; however, note that at 2: > 2 our model, built upon the 
SAM by G04, still performs quite well by construction. We 
will also pursue the analysis of galaxy statistics at multi- 
ple wavelengths; this should allow us to understand better 
the interplay between the processes involved within our sce- 
nario. Finally, we will discuss more extensively the structural 
properties of the discs emerging from our model, that will 
constitute testbeds for the next generation of SAMs. 
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